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Resumen 


Uno de los parámetros más difíciles de estimar al modelar las redes de distribución 
de agua potable es el consumo doméstico. Se ha demostrado que éste sigue un 
proceso estocástico y que es posible caracterizarlo a través de pulsos rectangulares 
con cierta intensidad, duración y apego a una frecuencia de arribo; todo ello a 
través de esquemas estocásticos, como el Modelo de Pulsos Rectangulares de 
Neyman-Scott (Neyman-Scott Rectangular Pulses Model, NSRPM; Neyman y Scott, 
1958). El esquema NSRPM se basa en la solución de un problema de optimización 
no lineal, que involucra momentos teóricos que representan las series sintéticas 
(equiprobables) y los momentos observados (mediciones de campo), que establece 
estadísticamente la serie medida. Se han publicado la metodología y validación 
inicial del proceso aplicado al consumo doméstico; sin embargo, no se ha llevado 
a cabo su validación en una red de distribución real y su comparación con los 
métodos tradicionales. En el presente artículo se comparan los resultados obtenidos, 
empleando simulaciones de series sintéticas con carácter estocástico, producto del 
esquema NSRPM aplicado a la determinación de presiones, caudales y fugas; con 
los obtenidos con el método tradicional que utiliza la Curva de Variación Horaria 
de la Demanda, y con mediciones de presión y caudal hechas en el sector Humaya, 
en Culiacán, Sinaloa, México. 
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EPANET* (Todini y Pilati, 1988), InfoWorks* 
y ScadRED*, entre otros. En otras palabras, lo 


Introducción 
En los últimos años se han desarrollado varios anterior se traduce en asignar el comporta- 


programas de cómputo, a fin de conocer en miento del consumo doméstico a través de una 


detalle el comportamiento hidráulico de los 
sistemas de distribución de agua potable. 

Una de las variables hidráulicas de estos 
modelos es el consumo doméstico, mismo que 
ha sido idealizado como una variación tempo- 
ral horaria. A esto se le conoce como “Curva de 
Variación Horaria de la Demanda (CVHD””, 
y es utilizada en programas comerciales de 
modelación de redes de agua potable, como 


curva de variación horaria. 

Sin embargo, la realidad de las tomas 
domiciliarias dista mucho de esta representa- 
ción, pues los consumos domésticos tienden a 
presentar un carácter estocástico (Buchberger ef 
al., 2003; Alvisi et al., 2003; Alcocer-Yamanaka, 
2007) y no permanente, sobre todo en escalas 
de tiempo del orden de segundos. Por ello, en 
los últimos años se han desarrollado técnicas 
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con un enfoque estocástico para representar 
el consumo doméstico, tales como el Método 
de Pulsos Rectangulares de Poisson (en inglés 
Poisson Rectangular Pulse Method, abreviado 
como PRP method) (Buchberger y Wu, 1995; 
Buchberger et al., 2003) y el Método de Pulsos 
Rectangulares de Neyman-Scott (en inglés 
Neyman-Scott Rectangular Pulse Method, 
abreviado como NSRPM) (Alvisi et al., 2003; 
Alcocer-Yamanaka et al., 2008a, 2008b). De 
estas técnicas, el método PRP maneja la 
demanda instantánea con intervalo de registro 
de un segundo; mientras que NSRPM puede 
trabajar con diferentes intervalos de registro, 
al considerar una desagregación temporal de 
la demanda. 

Estos métodos han tenido aplicación 
principalmente en el campo de la hidrología, 
en la generación de series sintéticas que 
representan lluvias o tormentas, según el 
intervalo y la duración proyectada. Las series 
generadas tienen parámetros estadísticos 
similares o idénticos a las series originales 
(observadas), tales como la media, varianza, 
covarianza y la distribución de probabilidad. 

La ventaja del NSRPM es que trabaja con 
diferentes intervalos de registro de los datos 
observados; representan la evolución de los 
consumos domésticos y tienen un mayor 
apego al comportamiento real de los usuarios 
del sistema. 

La estimación de los parámetros y, por 
ende, la generación de series sintéticas permite 
minimizar la información que sería necesario 
recabar con la medición de los consumos 
domiciliarios. 

En la bibliografía existente sólo se han 
registrado dos aplicaciones del esquema 
Neyman-Scott Rectangular Pulses Model (NSRPM) 
que intentan representar consumos domésticos 
(Alvisi et al., 2003; Alcocer-Yamanaka 2007). 


Sitio de aplicación 
Los esquemas determinístico y estocástico se 


aplicaron en el sector Humaya, en la ciudad 
de Culiacán, Sinaloa (ver figura 1), debido a 


la suficiente información de campo con que se 
cuenta (Alcocer-Yamanaka y Tzatchkov, 2002, 
2003, 2004; Alcocer-Yamanaka et al., 2004, 2007, 
2008a, 2008b; Tzatchkov et al., 2004, 2005), entre 
la que destaca lo siguiente: mediciones de 
presión y caudal en fuentes de abastecimiento 
y puntos internos de la red de distribución 
de agua; nivel en el tanque de regularización; 
calidad del agua en fuentes de abastecimiento 
y puntos dentro de la red de distribución (cloro 
residual, carbono orgánico total (COT), pH, 
temperatura); y consumo intradomiciliario en 
69 domicilios, con intervalo de registro de un 
minuto y duración promedio de tres días. 


Fuentes de abastecimiento 


La zona tiene dos captaciones empleadas 
como fuentes de abastecimiento: la primera se 
compone de un solo pozo, que aporta un gasto 
promedio de 511/s, y la segunda de una batería 
de ocho pozos, con capacidad máxima de 200 
litros por segundo. 


Estructuras de regularización y control 


Existen dos tanques de regularización: uno 
con capacidad de 3 000 m* y una cota de 82.63 
msnm, y otro con capacidad de 2 000 mY y cota 
de plantilla de 80.00 msnm. 


Población de la zona 


La población en el año 2005 era de 85 483 
habitantes; esto se determinó con base en el 
número de tomas (20 353 tomas en total) 
de cada colonia incluida dentro del sitio 
de aplicación del modelo, y el índice de 
hacinamiento por toma, que es de 4.20 hab / 
toma, de acuerdo con la información de la Jun- 
ta Municipal de Agua Potable y Alcantarillado 
de Culiacán (JAPAC). 


Fugas 


Según reportes del organismo operador, en la 
zona se tiene un porcentaje de pérdidas físicas 
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Figura 1. Representación del área de estudio, sector Humaya, en la ciudad de Culiacán, Sinaloa, México. 


del 30%, las cuales se originan, en su mayor 
parte, en tomas domiciliarias. La empresa JBS 
Associates Inc. realizó estudios de recupera- 
ción de agua durante los años 2000 y 2001 (BS, 
2000; JBS, 2001), que validan el porcentaje de 
fugas mencionado. 


Modelo utilizado 


Los datos geométricos de la red de abaste- 
cimiento de agua potable, que incluyen todas 
las tuberías (2” a 18” de diámetro) y otros datos 
necesarios para la modelación hidráulica, como 
la población y el consumo, se introdujeron en 


el programa de cómputo EPANET? (ver figura 
2). 

Una de las ventajas de este programa es 
su capacidad y versatilidad para analizar 
periodos extendidos, dado que en el esquema 
estocástico se modeló el lapso de una semana, 
es decir, 168 horas. 


Sitios seleccionados para comparar 
resultados de los dos esquemas con las 
mediciones en la red 


Los resultados obtenidos con los esquemas 
determinístico y estocástico se compararon con 
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Figura 2. Red de distribución de la zona Humaya digitalizada dentro del sistema EPANET". 


las mediciones de campo, de esta forma fue 
posible establecer las ventajas y desventajas de 
ambos. En la figura 3 se muestra la ubicación 
de los nodos y tramos analizados dentro del 
sector Humaya, Culiacán, Sinaloa, aunque por 
razones de extensión de este trabajo, sólo se 
discutirán algunos. 


Enfoque determinístico. Curva 
tradicional de la demanda 


La curva de variación horaria es una ideali- 
zación del consumo obtenida a partir de 
mediciones domésticas y comerciales en 
sectores de la red aislados hidráulicamente, 
llamados distritos hidrométricos. La Comi- 
sión Nacional del Agua la obtuvo mediante 
el registro de los consumos medidos 
aleatoriamente en los distritos con tres 
diferentes niveles socioeconómicos y uno 
comercial, en varias ciudades, comparándolos 
con el caudal medido a la entrada de los 


distritos (ver figura 4). 
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Figura 3. Ubicación en la red de nodos y tramos analizados. 


Es importante mencionar que la metodo- 
logía empleada para el trazo de la curva 
incluye los consumos de los usuarios y las 
fugas dentro de la red analizada. Asimismo, se 
puede observar que la curva está suavizada, 
lo cual, como se mostrará más adelante, 
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Figura 4. Curva de Variación Horaria de la Demanda para la república mexicana (Tzatchkov, 2007). 


no ocurre en la realidad, al compararla con 
mediciones continuas de los consumos de una 
casa-habitación. 


Presiones y caudales calculados 


Se presentan las presiones calculadas con el 
programa EPANET, en el nodo 165 (figura 5) 
(ver ubicación en la figura 3). 

Se puede observar que la variación de 
presión tiene un comportamiento cíclico, 
dado que se mantiene la demanda base y el 
patrón de comportamiento es introducido por 
medio de la Curva de Variación Horaria de la 
Demanda. Las presiones máximas y mínimas 
medidas en metros de columna de agua (mca) 
son 26.00 y 22.00 mca, respectivamente. 

La figura 6 muestra el área de influencia del 
tramo 2957, de 12 pulgadas de diámetro, que 
suministra una zona dentro del sitio analizado; 
la figura 7 ilustra el comportamiento de los 
caudales obtenidos con el mismo programa. En 
éste se tiene también un comportamiento cíclico 
y caudales máximos de 28.0 1/s y mínimos de 
10.0 1/s (aproximadamente), todos positivos. 


Conclusiones de este esquema 


Se observa que las curvas de comportamiento 
de caudales y de presiones son cíclicas y sua- 
vizadas debido a la escala horaria que establece 
la Curva de Variación de la Demanda. 

En el tramo analizado se observó que el 
caudal tiene signos siempre positivos; el hecho 
de que no haya cambio de signo en el caudal 
significa que el flujo sigue el mismo trayecto de 
“entrada” y “salida” del tramo. 


Enfoque estocástico 


El reconocimiento de que los consumos se 
producen aleatoriamente ha motivado a al- 
gunos investigadores (Kiya y Murakawa, 
1989; Buchberger y Wu, 1995) a proponer 
que la ocurrencia de consumos a lo largo del 
tiempo sigue un proceso de Poisson. Esta 
hipótesis fue verificada por Buchberger et al. 
(2003). Dicho proceso resulta no homogéneo 
debido a que los consumos a lo largo del día 
varían considerablemente de unas horas a 
otras. Cada consumo individual de agua se 
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Figura 5. Presión obtenida en el nodo 165 para el escenario, empleando la CVHD y el programa EPANET. 
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Figura 6. Área de influencia del tramo 2957 y su salida a través del tramo 22. 


representa como un pulso rectangular con una 
altura (intensidad) x y un ancho t (duración) 
aleatorios. 
Recientemente se han desarrollado 
modelos de simulación de consumo como el 


Poisson Rectangular Pulses (PRP) (Buchberger 


et al., 2003), que permiten generar series de los 
consumos bajo ciertos criterios estocásticos, 
que se basan en los siguientes parámetros 
básicos: tasa de llegada o frecuencia de la 
ocurrencia de los pulsos individuales (WM), 
intensidad promedio de los pulsos (Hu), 
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Figura 7. Caudales calculados para el tramo 2957, con la CVHD y EPANET. 


varianza de esa intensidad Var(u ), duración 
promedio de los pulsos (n) y la varianza de 
esa duración Var(n). La obtención de estos 
parámetros se ha realizado generalmente 
con base en mediciones de consumo con 
intervalo de registro de un segundo. Aunque 
realizar mediciones con intervalo de registro 
de un segundo tiene la ventaja de obtener 
directamente la evolución real del consumo 
doméstico, se requiere de sofisticados equipos 
de medición y almacenamiento de datos, y 
de un elevado esfuerzo computacional en el 
análisis de los datos generados (Buchberger 
et al., 2003). Estas razones limitan métodos 
como el PRP y dificultan su aplicabilidad en 
proyectos de investigación en la actualidad, 
así como el desarrollo de proyectos de 
diseño y análisis de las redes de distribución 
en el futuro, en especial en países en vías 
información e 
infraestructura generalmente es limitada. 


de desarrollo, donde la 

En los últimos años se han desarrollado 
técnicas orientadas a la estimación indirecta 
de los parámetros 1, u,, Var(u ), n y Var(n), de 
registros de la demanda con intervalos más 
largos, destacando la desagregación espacial 
y temporal (Alcocer-Yamanaka et al., 2006, 
2007, 2008a, 2008b, 2009; Guercio et al., 2001; 


Rodríguez-Iturbe et al., 1984, 1987). En esas 
técnicas, la estimación de los parámetros 
se basa en el planteamiento de una función 
objetivo que expresa la relación entre los 
momentos estadísticos de la serie observada 
y los momentos teóricos del modelo de 
consumo, que se minimiza a través de técnicas 
de programación no lineal, obteniéndose 
a partir de ello los parámetros deseados. 
Nadimpalli y Buchberger (2003) realizaron 
una comparación entre esas técnicas aplicadas 
al problema de estimación de los parámetros 
con base en ejemplos. En todas ellas se supone 
conocida la variación de la demanda en una 
tubería que abastece cierto número de casas. 
Las técnicas se diferencian entre sí por el 
tipo de distribución de probabilidad que se 
asume para gobernar el comportamiento de 
algunos parámetros, por ejemplo, la duración 
e intensidad promedio de los pulsos; la tasa 
de llegada de los pulsos, y por el proceso 
estocástico que se emplea como base para 
formular los momentos teóricos involucrados 
(Rodríguez-Iturbe et al., 1984). 

El esquema de Neyman-Scott ha sido 
empleado previamente en la generación del 
consumo por Alvisi et al. (2003), y Alcocer- 
Yamanaka et al. (2007, 2008a, 2008b, 2009a, 
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2009b). Los momentos de segundo orden 
del proceso agregado Y(” son los siguientes 
(Entekhabi y Bras, 1990): 


E[x/0] — Hch (1) 


Var [y] = A (nh-1+e "| CS 


+E|c?-C] ¿a (2) 


n? 
2 
nespe] lo qui | (3) 
al A 
All-e py He Clu%, Am BUD > 1 


En donde 1* representa el tiempo promedio 
entre dos eventos; [*, el tiempo promedio 
entre cada pulso individual y el origen del 
evento; 1], la duración promedio de los pulsos; 
u, la intensidad promedio de los pulsos, y h 
es el intervalo de agregación / desagregación 
analizado. 

Definidas las expresiones del esquema de 
NSRPM se formula la función objetivo: 


De donde F'1, F2,... F'n son los valores de 
los momentos observados; es decir, la media, 


varianza y correlación lag-1, entre otros. Por 
su parte, Fl, F2, F3,...Fn son los momentos 
teóricos, funciones del vector de parámetros, 
E (A, HL, 1, n, B). Para este caso de aplicación en 
consumos domésticos se consideró n = 3, que 
representa la media, varianza y covarianza en 
la ecuación (4). 

En la formulación del esquema de Ney- 
man-Scott se debe establecer el intervalo de 
análisis buscado para la ejecución del esquema 
de optimización (para este caso en particular, 
el tiempo será de un minuto). Enseguida se 
realiza, a través de técnicas de programación 
matemática no lineal (NLP) (método de 
gradiente conjugado con derivadas centrales 
y estimación cuadrática), la minimización de 
la función objetivo, con lo cual se obtienen 
valores para cada uno de los parámetros del 
modelo (llamadas variables de decisión en 
un esquema de optimización). 

Para aplicar el esquema de Neyman-Scott 
en un modelo de red de agua potable es nece- 
saria la asignación de patrones estocásticos 
generados de manera aleatoria a la demanda 
en cada nodo del modelo (cada nodo contiene 
cierta cantidad variable de viviendas), en 
función del nivel socioeconómico, lo que 
implica introducir patrones de consumo con 
una longitud de 10 080 datos, correspondien- 
tes a cada minuto transcurrido durante 
una semana. La asignación de los patrones 
estocásticos debe considerar también el nivel 
socioeconómico de las viviendas, para lo cual 
se dividieron en tres grupos (bajo, medio y 
alto), y se generaron patrones estocásticos 
diferentes para cada grupo. 

Para tomar en cuenta el hecho de que la 
demanda de agua potable es diferente durante 
el día, las series generadas, con sus parámetros 
estadísticos, se dividieron en tres bloques 
horarios, como se describe más adelante en 
este artículo. Inicialmente se determinaron 
los parámetros necesarios para la generación 
de las series sintéticas en las 69 casas- 
habitación en que se registró la variación de 
la demanda temporal, después se generaron 
cincuenta series sintéticas, con la finalidad de 
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ensamblarlas para su validación. Entiéndase 
por “ensamble” al proceso de generar cincuenta 
series y posteriormente promediar los valores 
de los momentos estadísticos de todas éstas. 

Así, comparando el momento observado 
con el momento de las series ensambladas 
de cada bloque horario, se determinó si se 
aceptaban o no las series sintéticas corres- 
pondientes para ser utilizadas en el esquema 
de simulación estocástico. En otras palabras, 
este proceso establece que en los casos donde 
los valores entre los momentos (observados 
y ensamblados) difieran, las series sintéticas 
generadas serán consideradas como no válidas 
para el patrón analizado y serán desechadas; 
en su lugar se generarán nuevas series. 
Caso contrario, cuando la diferencia entre 
los momentos tienda a cero, el proceso y las 
series sintéticas generadas son considerados 
válidos. 

Posteriormente se generaron 69 curvas de 
patrones de consumo sintéticos, asignando 
aleatoriamente bloques horarios de la serie 
generada para formar un patrón sintético de 
una semana de longitud de lunes a domingo, 
distribuyendo aleatoriamente los 69 patrones 
sintéticos a cada casa-habitación, en función 
de su nivel socioeconómico. 


Tratamiento del registro de datos de 
consumos domésticos 


Para el tratamiento de los datos registrados 
se realizó un promedio de los consumos, 


separando los registros en bloques horarios 
de los diferentes días de la semana, en función 
del comportamiento típico de la demanda de 
agua potable (ver cuadro 1), donde puede 
observarse la separación de los días labora- 
bles, fines de semana y días festivos. La 
separación en bloques de los datos corresponde 
a la disparidad numérica que presentan los 
valores de los momentos observados en los 
distintos horarios. Ello facilita la obtención 
de una solución óptima durante el proceso 
de optimización y, por ende, la generación 
de series sintéticas con mayor apego a la 
realidad. 

Una vez separados los datos registrados en 
campo en los bloques marcados en el cuadro 
1 se obtienen tres bloques y de éstos, ocho 
subbloques horarios: cuatro para los días de 
lunes a viernes, dos para el sábado y dos para 
el domingo. 

Una vez obtenidas las series promediadas 
para cada bloque horario se determinan sus 
momentos, de media, varianza, covarianza y 
volumen acumulado. Posteriormente, para 
obtener los parámetros estadísticos adecuados 
y válidos dentro del esquema Neyman-Scott 
es necesario acotar el espacio de solución del 
esquema de optimización, empleando rangos 
de búsqueda a partir de la reducción de los 
intervalos en el orden de magnitud de las 
variables de decisión. Finalmente se generan 
las series sintéticas que serán implementadas 
dentro del programa de dominio común 
EPANET. 


Cuadro 1. Bloques horarios para el análisis de consumos registrados en función del comportamiento 


típico de la demanda. 


Del primer domingo de abril al sábado anterior al último domingo de octubre 


Día de la semana 1* bloque 2" bloque 3* bloque 
6:00 - 20:00 
Lunes a viernes 0:00 - 6:00 20:00 - 22:00 
22:00 - 24:00 
Sábado 0:00 - 7:00 7:00 - 24:00 
Domingo y festivo 0:00 - 19:00 19:00 - 24:00 
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Generación de series sintéticas y 
determinación de rangos de búsqueda en 
el esquema de optimización para obtener 
los parámetros del modelo Neyman- 
Scott 


Una vez obtenidos los momentos observados 
para los domicilios donde se hicieron las medi- 
ciones, mediante el esquema de optimización 
se determinan los parámetros estadísticos (A, 
x, UM, c, y, BP) que intervienen en los momentos 
teóricos representados en la ecuación (4) 
(función objetivo), para su posterior intro- 
ducción en el modelo NSRPM. 

La generación de las series se hizo con base 
en el modelo de dominio público contenido 
en el Rainfall Data Modelling Portal (RDMP) 
(Mellor, 2007). Dado que se trata de un evento 
estocástico, es importante señalar que se 
deberá realizar cierto número de simulaciones 
dentro del NSRPM, considerando una semilla 
de generación de números aleatorios diferente 
entre una simulación y otra. Finalmente, se 
comparan las series sintéticas obtenidas a 
partir del NSRPM con la serie original medida 
en campo. 

Los rangos de búsqueda inicialmente utili- 
zados para el cálculo de la función objetivo y 
de los parámetros para cada vivienda fueron 
los indicados en el cuadro 2. 

Una segunda aproximación de los rangos 
de búsqueda de los parámetros mencionados 
se hizo usaron datos reportados por Feliciano 
(2005), se usaron valores de frecuencia A” 
desde 0.0689 min? (14.51 minutos) hasta 
los 0.04305 min” (23.23 minutos). Con éstos 
se determinó un rango de búsqueda de 1 
min”? (1 minuto) hasta 0.0404 min” (24.75 
minutos), para tener un margen más amplio de 
localización en el esquema de optimización, 
pues se observó que al dejar un rango menor 
O igual a los 24.75 minutos, ocasionaba que el 
resultado en la función objetivo se alejara de 
cero. Es importante mencionar que los rangos 
reportados por Feliciano (2005) corresponden a 
una zona aledaña a la actual y el intervalo de 
registro empleado fue de un segundo. 


Cuadro 2. Rango de búsqueda inicial para los parámetros 
del modelo Neyman-Scott. 


Parámetro Desde Hasta 
A*(min”) 1 0.016667 
B” (min”) 1 0.05 
n* (min”) 10 0.007 
u, (1/min) 1 8 

E(c) 1 20 


El siguiente parámetro que se ajustó fue 
la intensidad promedio de los pulsos «4. Los 
valores medios de la intensidad de los pulsos 
considerados fue de 1 a 8 1/min, sin embargo 
tales valores “forzaban” en cierta forma al 
esquema de optimización, teniéndose como 
resultado valores de volumen consumido 
muy altos (irreales). Después de varias pruebas, 
se determinó que había que reducir el valor de 
8 a 6 1/min, para hacer menor la intensidad 
promedio y con esto disminuir la intensidad 
de los pulsos en la serie sintética, lo cual 
resultó en que también la varianza observada 
se aproximase de mejor manera a la varianza 
de las series sintéticas (ver cuadro 3). 

En el cuadro 4 se muestran los resultados 
únicamente de tres de las viviendas (el 
desglose se omite por motivos de extensión del 
manuscrito). 

Como se mencionó anteriormente, con 
los parámetros calculados se procedió a la 
generación de las cincuenta series sintéticas con 
datos a cada minuto, por bloque horario de 
cada vivienda, generando aproximadamente 27 
000 series sintéticas, que cubren los consumos 


Cuadro 3. Rango de búsqueda definitivo para obtener los 
parámetros del modelo Neyman-Scott. 


Parámetro Desde Hasta 
A*(min”) 1 0.0404 
B” (min”) 1 0.05 
n” (min”) 10 0.007 
u, (1/min) 1 6 

E(c) 1 20 
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de una semana para las 69 viviendas analizadas. 
Los resultados obtenidos en el cálculo del 
esquema de optimización para la generación 
de los parámetros que requiere el modelo 
Neyman-Scott se presentan en Hernández 
(2009). 

En cada nodo del modelo existe una 
cantidad no uniforme de casas asignadas, 
por lo que para cada casa se asigna un patrón 
estocástico y una demanda media con base en 
éste. Se consideró una duración de los siete 
días de la semana, y se asignaron las demandas 
en los nodos en función de la cantidad de 
predios que cubría cada uno de ellos. Las 
demandas asignadas se obtuvieron de los 69 
patrones sintéticos de consumo generados 
previamente. Es importante mencionar que los 
patrones sintéticos corresponden al consumo 
de las casas-habitación, de tal forma que en 
EPANETO se introduce cada patrón sintético 
a la casa correspondiente. 

De manera aleatoria se asignaron los 
69 patrones sintéticos en función del nivel 
socioeconómico de las casas-habitación. Los 
patrones contienen 10 080 datos, que corres- 
ponden a los consumos con intervalo de 
registro de un minuto y duración de siete 
días. 

Enseguida se muestran los resultados para 
este esquema en los mismos nodos y tramos 
que se consideraron en el determinístico, y 
con ello se realiza una comparación de los 
resultados de ellos y con mediciones de campo. 
Vale la pena señalar que el caudal asignado 
para generar en el programa EPANETO la 
curva característica sintética del bombeo en 
la batería de pozos fue de 104.22 litros por 
segundo. 


Comparación de resultados entre los 
modelos estocástico y determinístico 


Presiones y caudales calculados 
Se presentan las presiones calculadas con el 


programa EPANET en el nodo 165 (figura 
8) (ver ubicación en la figura 3). Se observa 


una fluctuación diaria debido a la variación 
estocástica del patrón de las series sintéticas, 
iniciando el día domingo. Las presiones 
máximas y mínimas medidas en metros de 
columna de agua (mca) son 26.34 y 2.69 mca, 
respectivamente. 

Continuando el análisis se revisó el 
comportamiento del flujo en el mismo tramo 
examinado con el esquema determinístico 
(tramo 2957). En la figura 9 se muestra el flujo 
en el tramo (suministro a una zona del sitio 
analizado) de 12 pulgadas de diámetro, donde 
se obtuvieron caudales en un solo sentido 
de la tubería, con valores entre 106.08 1/s y 
cercanos a cero litros por segundo. 

El comportamiento de los caudales y 
presiones es muy variable bajo este esquema. 
Los cambios repentinos son ocasionados 
por la aleatoriedad con que se formaron los 
patrones de consumo, lo que provoca que 
en ciertos instantes (del orden de minutos) 
existan consumos altos y en el siguiente 
la demanda disminuya hasta el punto de 
llegar a ser nula. Todas estas variaciones 
las absorben el tanque de regulación y los 
equipos de bombeo. 

Los caudales que se presentan en los 
nodos en este esquema pueden ser positivos o 
negativos a lo largo del tiempo, lo que indica 
cambios de dirección del flujo; la ausencia 
de flujo, por su parte, representa tiempos 
de estancamiento O tiempos de residencia 
elevados del agua dentro de las tuberías. 

Es importante señalar que los patrones 
estocásticos se pueden variar utilizando un 
bloque cualesquiera de las series sintéticas 
generadas; en este caso se utilizó aleatoria- 
mente el bloque 33. 

En el cuadro 5 se muestra el resumen 
de resultados de presión para el nodo 165, 
considerando ambos enfoques. 

Las presiones máximas y mínimas que se 
presentan en el escenario estocástico ocurren 
en lapsos hasta de un minuto, siendo éste 
el tiempo de análisis que se introdujo en la 
simulación. Vale la pena explicar que las 
presiones mínimas que se presentan en los 
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Figura 8. Presión obtenida en el nodo 165 para el escenario estocástico. 
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Figura 9. Caudales calculados para el tramo 2957 en el escenario estocástico. 


resultados de la simulación con el modelo 
Cuadro 5. Resumen de resultados de presión en ¿os ¿ ES ; 
. estocástico, en particular la presión negativa 

el nodo 165 para ambos escenarios. 


que se obtiene, corresponden a las premisas 


NT an en que se basa el programa EPANET, es decir, 

a Presión máxima Presión O . . 
Escenario | Nodo a ia (En un análisis de periodos extendidos, donde 
VE 166 dE 4158 las variaciones de gasto y presión son lentas, 
l - pero pueden ser no realistas en este caso de 
pa o ME ES pt demanda estocástica bruscamente variable. Ello 


indica que la modelación debería hacerse con 
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un modelo dinámico más refinado, capaz de 
representar variaciones bruscas de las variables 
hidráulicas, pero la discusión de tal modelo 
está fuera del alcance del presente artículo. 

En relación con los flujos, en el cuadro 
6 se observa que en el esquema estocástico 
existen caudales mayores a los obtenidos con el 
esquema utilizando la CVHD. 

Finalmente se observa que las presiones 
y los caudales resultan diferentes con ambos 
esquemas debido a la observación hecha ante- 
riormente. 


Comparación de presiones y caudales 
obtenidos con los esquemas planteados y 
las mediciones de campo 


En campo se efectuaron mediciones de presión 
y caudal en varios nodos y tramos en los cuales 
se aplicaron los esquemas determinístico y 
estocástico; aquí se presentarán sólo algunos 
por razones de espacio. En la figura 10 se 
muestra la comparación entre las presiones 
medidas en el punto ubicado en la colonia 
“Santa Fe” y las obtenidas con ambos esquemas 
en el nodo 165. Se observa que las mediciones 
y el esquema estocástico representan mejor el 
comportamiento del consumo por pulsos, y 
que si se obtienen curvas promedio en ambos 
casos, la coincidencia es mayor que con la 
CVDH. 

Como en el caso de las presiones, las 
mediciones de caudal y el esquema estocástico 
reflejan mejor el patrón de consumo por pulsos 
(ver figura 11). 

Otra ventaja del enfoque estocástico es 
que permite estimar las fugas en la red. Por 
la forma en que fue obtenida la CVHD, ésta 
incluye las pérdidas físicas; si se compara con 


las mediciones, puede notarse que la CVHD 
se ubica por encima de la curva que representa 
la demanda de los usuarios. 


Conclusiones 


En este artículo se muestra la aplicación de 
una metodología que considera conceptos 
estocásticos para determinar patrones de 
consumo doméstico. Se emplea el esquema 
de Neyman-Scott, que al introducirse en 
un modelo de simulación hidráulica per- 
mite obtener resultados más parecidos al 
funcionamiento de una red de agua potable. 

Actualmente el enfoque empleando 
la CVHD resulta un apoyo sencillo para 
los ingenieros e investigadores dedicados 
al análisis y el diseño de los sistemas de 
distribución de agua potable. Sin embargo, esta 
curva tiende generalmente a sobredimensio- 
nar y, en algunas ocasiones, subdimensionar 
la infraestructura hidráulica, producto de los 
resultados obtenidos en las simulaciones y, por 
ende, diseñadas de esta forma. Se recomienda 
al proyectista tener en consideración esta 
observación y si de todas formas usa la CVHD, 
obtenerla por mediciones en campo para cada 
sitio analizado. 

El resultado de este trabajo es el inicio en la 
aplicación de este enfoque tendiente a conver- 
tirse práctico y sencillo para los ingenieros 
e investigadores dedicados a los sistemas 
de distribución de agua potable. Ello podría 
concretarse a través de la incorporación de 
esta metodología a través de un módulo dentro 
de programas comerciales como EPANET. 

Finalmente, en el ámbito local, la imple- 
mentación del modelo desimulación hidráulica, 
considerando el consumo apoyado en un 


Cuadro 6. Resumen de resultados de caudales en tramos para los dos escenarios. 


Escenario Tramo Colonia ida Q máx (l/s) Observaciones 
CVHD 2957 12 28.13 Sin cambio en sentido del flujo 
Infonavit Humaya 
Estocástico 2957 12 106.08 Sin cambio en sentido del flujo 
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Figura 10. Comparación de las presiones registradas en el nodo 165 con los esquemas utilizados. 
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Figura 11. Comparación con caudales medido, resultante con CVHD y estocástico, en el tramo 2957. 
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enfoque estocástico, permitirá al organismo 
operador emplearlo como toma de decisiones 
ante posibles fuentes de abastecimiento exter- 
nas, dado que es una zona de incremento 
poblacional elevado, así como en el diseño de 
un programa de rehabilitación y sustitución de 
tuberías debido a los problemas de incrustación 
de fierro y manganeso que sufre el área de 
estudio. 

Se recomienda como trabajo futuro auto- 
matizar el proceso de generación de las series 
de demanda estocástica, estableciendo un 
esquema de Monte Carlo en el proceso de 
simulación. 
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Abstract 


ARREGUÍN, El., ALCOCER-YAMANAKA, V.H. € HERNÁNDEZ-PADRÓN, D.S. 
Modeling of drinking water distribution networks using deterministic and stochastic 
approaches. Water Technology and Sciences, formerly Hydraulic engineering in Mexico 
(in Spanish). Vol. 1, No. 4, October-December, 2010, pp. 119-136. 


Home water consumption is one of the most difficult parameters to estimate while modeling 
drinking water distribution networks. It has been shown that this parameter follows a 
stochastic process, and that it is possible to characterize it with rectangular pulses —having a 
certain intensity, duration, and adherence to an arrival frequency —using stochastic models, 
such as the Neyman-Scott Rectangular Pulses Model (NSRPM). This model is based on the 
solution of a nonlinear optimization problem involving theoretical moments that represent 
the synthetic (equiprobable) series, and observed moments (in-field measurements) that 
statistically establish the measured series. The methodology and the initial validation of the 
process applied to home consumption have been published; however its validation in an actual 
distribution network has not been made, nor has it been compared with conventional methods. 
This paper compares the results obtained by using simulations of synthetic series which are 
stochastic in nature, product of the NSRPM model applied to the determination of pressures, 
flows, and leaks, with those obtained with the traditional method that uses demand time- 
variation curves and with pressure and flow measurements made in the Humaya sector, in 
Culiacán, Sinaloa, Mexico. 


Keywords: modeling, water distribution networks, demand. 
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